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We study inclusive di-pion decays using a sample of 108 x 10*' T(3S') events recorded with the 
BABAR detector. We search for the decay mode T(3S') it'^tt~ hb{lP) and find no evidence for the 
bottomonium spin-singlet state hi,{lP) in the invariant mass distribution recoiling against the -K^n~ 
system. Assuming the hb{lP) mass to be 9.900 GeV/c^, we measure the upper limit on the branching 
fraction B[T{3S) ■k'^-k~ hb{lP)] < 1.2 x 10^**, at 90% confidence level. We also investigate the 
Xbj{2P) 7r+7r-X6j(lP), T(35) ^ 7r+7r"r(2S'), and T{2S) n+TT-T{lS) di-pion transitions 
and present an improved measurement of the branching fraction of the T{SS) — !> 7r^7r~T(25') decay 
and of the T(3S) — T{2S) mass difference. 



PACS numbers: 13.20.Gd, 14.40.Pq, 14.65. Fy 

Studies of bb (bottomonium) and cc (charmonium) 
bound states provide insight about inter-quark forces. 
The measurement of the hyperfine mass splitting be- 
tween triplet and singlet states in quarkonium systems 
discriminates between various models and tests lattice 
QCD and potential nonrelativistic QCD calculations fll]. 
Observation of the P-wave singlet ground state of char- 
monium, hc{lP), was recently confirmed and its mass 
precisely measured, yielding the hyperfine splitting for 
the charmonium IP states AMhf{lP)cc = {M{^Pj))cq- 
M{^Pi)cs = +0.08±0.18(stat.)±0.12(syst.) MeV/c^ 0, 
where (Af(^Pj)) is the spin-weighted average mass of 
the J = 0,1,2 ground states. The hyperfine split- 
ting AMhfilP),-, = (Mi^Pj)),-, - M{^P,\j, for bot- 
tomonium states is expected to be no more than 
a few MeV/c^ The ^Pj bb ground states are 



well-established, and their spin-weighted mass aver- 
age is (M(3Pj),5) = [M{xm{lP)) + 3M(xw(lP)) + 
5M(xfc2(l^'))]/9 = 9.89987 ± 0.00027 GeV/c2 [i]. The 
hb(lP), hereafter referred to as the h^, is expected to de- 
cay predominantly to ggg (57% branching fraction), 7776 
(41%), and 537 (2%), and its width is predicted to be of 
order 0.1 MeV [5^. 

We report, herein, a search for the /ib through 
the hadronic transition TiiS) — > 7r+7r^ft,b(lP). The 
CLEO Collaboration searched for the hi, in the reac- 
tions T[?>S) — > 7r"/if, and T(3S') — > n^n^hb, setting up- 
per limits at 90% confidence level (CL) for the branch- 
ing fractions B[r{iS) TT°hb] < 2.7 x lO'^ and Br = 
B\T{iS) 7r+7r~ft,;,] < 1.8 x 10~^, assuming the hb mass 
m{hb) to be 9.900 GeV/c^ The BaBar Collaboration 
recently reported evidence for the hb in T(3S') — >■ ir^hb. 
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7^67 decays Preliminary results of a search for 



the hh in the reaction e+e 



TT /if,, reporting the ob- 



servation of the hb meson, have been announced by the 
Belle Collaboration Theoretical predictions for Br 
span one order of magnitude. References 0-ll| predict 
a branching fraction between 2.2 x 10^^ and 8.0 x 10~^, 
while Ref. 12| predicts a rate of lO"** or smaller. 

The data sample used in this study was collected with 
the BaBAR detector at the PEP-II asymmetric-energy 
e+e" storage rings at the SLAC National Accelerator 
Laboratory. It consists of 25.6 fb~^ of integrated lumi- 
nosity collected at a e"^e~ center-of-mass (cm.) energy of 
10.355 GeV, corresponding to the mass of the T{3S) reso- 
nance. The number of recorded T{3S) events is 108 x 10^. 
An additional sample of 2.5 fb~^ recorded at the T{3S) 
energy ("10%" sample) and a 2.6 fb""'^ sample collected 
30 MeV below the T{3S) resonance ("off-peak" sample) 
are used for background and calibration studies. 

The momenta of charged particles are reconstructed 
using a combination of a five-layer double-sided silicon- 
strip detector and a 40-layer drift chamber, both oper- 
ating in the 1.5-T magnetic field of a superconducting 
solenoid. Photons are detected using a CsI(Tl) electro- 
magnetic calorimeter, which is also inside the magnet 
coil. Charged hadron identification is achieved through 
measurements of particle energy loss in the tracking sys- 
tem and the Cherenkov angle obtained from a detector 
of internally reflected Cherenkov light. 

The 7r"'"7r~ pairs are selected from oppositely-charged 
tracks that originate from the e~^e^ interaction region in 
hadronic events, hence excluding tracks arising from a 
photon conversion or the decay of a long-lived particle. 
We search for an hj, signal using a fit to the spectrum of 
the mass rrifj recoiling against the tt'^7t~ system, defined 
by: 



(M: 



T(3S) 



(1) 



where E*^ and P*^ are, respectively, the measured tttt 
energy and momentum in the cm. frame. 

The hh signal is expected to appear as a peak in 
the mji distribution on top of a smooth non-peaking 
background from continuum events {e~^e~ — qq with 
q = u,d, s, c) and bottomonium decays. Several other 
processes produce peaks in the recoil mass spectrum 
close to the signal region. Hadronic transitions T{iS) 
7r"'"7r~T(2S') (hereafter denoted T'^"*'^) produce a peak 
centered at the Y{2S) mass to[T(2S')] = 10.02326 ± 
0.00031 GeV/c^ The cascade process T(3S') 

T{2S)X, T{2S) -> 7r+7r-T(lS') (T^^^) results in a peak 
centered at 9.791 GeV/c^. The peak is offset from the 
Y{1S) mass by approximately the T{3S) to Y{2S) mass 
difference. Doppler shift and broadening further affect 
the position and width of this peak. When the T{2S) 
parent in decays is produced through the T^^"^ 

channel, a pion from the T{3S) decay can be combined 



with an oppositely-charged track from the T{2S) de- 
cay to produce a broad distribution centered around 
9.9 GeV/c^. The r(2S) is also produced through the 
initial-state radiation (ISR) process e+e" jjsrT{2S) 
C^isr)- Of the nine possible T{3S) XbJ'i^Ph, 
Xb.]'{2P) TT+n-XbjilP) decay chains (xf^^) M, only 
those for J' = J = {1,2} have been reported [J, |l5(; 
these should generate two narrowly separated peaks near 
9.993 GeV/c^, while the contributions with J' ^ J or 
with ,7 = are expected to be negligible. 

Selection criteria are chosen by optimizing the ratio 
S/^/B between the expected hh signal yield (5) and the 
background [B). The signal sample for the optimization 
is provided by a detailed Monte Carlo (MC) simulation 
based on Geant4 (16|], EvtCen pj], and Jetset [H], 
while the background sample is obtained from the 10% 
sample, which is not used for the extraction of the signal. 
The natural width of the hh meson, which is predicted to 
be negligible in comparison with the experimental reso- 
lution in rriB, (0.009 GeV/c-^ r.m.s.), is set to zero in the 
simulation. 

Since decays of the hh via three gluons or to rjh^ are ex- 
pected to exhibit a high track multiplicity, we require an 
event to have between 6 and 16 charged tracks, where the 
upper restriction reduces contributions due to random 
combinations of tracks. We further require the ratio of 
the second to zeroth Fox- Wolfram moment [l3| calculated 
using all charged tracks and unmatched neutral showers 
in the event to be less than 0.55. The total event energy 
in the laboratory frame must lie between 6 and 18 GeV, 
where the lower restriction reduces QED background. 

Events must contain two oppositely-charged tracks, 
each of which is identified as a pion. The pion identifi- 
cation efficiency depends on momentum and polar angle, 
and is typically about 98%. This requirement provides a 
rejection factor of order 50 against electrons. The vertex 
of each reconstructed pion pair must lie within 0.41cm 
and be less than 4(Tl from the nominal interaction point 
in the transverse plane, where is the uncertainty eval- 
uated on a candidate-by-candidate basis for the trans- 
verse fiight length L. We demand the probability for 
the vertex fit to be greater than 0.001. 

The phase-space distribution of decays extends up 
to mfl values of approximately 9.86 GeV/c^ and then 
rapidly decreases. To further suppress the background 
due to decays, we reject pairs of pious if their vertex 
is displaced from the nominal interaction point by more 
than 0.05 cm and 2(Tl in the transverse plane and if they 
satisfy cos a > 0.95, where a is the angle between the di- 
rection of the di-pion candidate momentum and its flight 
direction in the transverse plane. Candidates removed 
from the nominal sample that satisfy all other selection 
criteria constitute a i^Tg-enriched control sample. 

The selected data sample consists of approximately 
137 X 10^ 7r+7r~ candidates in the range 9.750 < mn < 
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10.040 GeV/c^, corresponding to an average of 2.4 se- 
lected di-pion candidates per event. The fit validation 
studies described below account for the effect of candi- 
date multiplicity. We evaluate the di-pion reconstruc- 
tion efficiency with MC events, by matching the recon- 
structed TT+TT^ pairs to the simulated pairs emitted in 
the bottomonium transition under study on an event-by- 
event basis. The hi, signal efficiency is e — 41.8% for 
m{hb) = 9.900 GeV/c2, with a [+6,-3]% variation of e 
over the m{hb) range [9.880,9.920] GeV/c^. A lower re- 
construction efhciency of 25.0% (16.7%) is found for the 



softer 



/r^ pairs produced in Xb ' (T^^^) transitions. 
For the T^^^ transition, an efficiency of 47.2% is ob- 
tained by averaging over the contributions from X = 77, 
tt'^tt'', and tt+tt". 

We perform a fit to the ma spectrum in the range 
9.750 < mR < 10.040 GeV/c^ with a model compris- 
ing eight components: non-peaking background, T^"*"^, 
Tf^jl, X6'\ TT+TT", and the K sig- 

nal. The mji distributions of the signal and background 
are parameterized using probability density functions 
(PDFs). We define a two-sided Crystal Bah (TCB) func- 
tion, which is a Gaussian for —ul < (x — xo)/a < an, 
and transitions to the power-law tail function f{x) [20| : 



\x - xo\ ^ nj 
a aj 



, (2) 



where and a are the mean and width of the Gaussian, 
and the subscript i = L {i = R) applies to values x < 
xq {x > xq). We model the signal component with a 
symmetric (q!l = a^, ul = nn) TCB shape. 

The and T^^^ peaks are described by sums of 

an asymmetric TCB shape and an asymmetric Gaus- 
sian. Contributions to T'^^^ from X = {tt+tt", tt^tt", 77} 
are modeled separately because of the different Doppler 
broadening. Their relative fractions and relative peak 
positions are fixed according to the world-average val- 
ues Q and the MC-simulated m_R spectrum, respectively. 
For each peak, the ratios of the widths of the TCB and 
Gaussian functions are fixed to the values found from fit- 
ting the corresponding MC spectrum. The PDF of the 
peaking background from ISR T{2S) production is pa- 
rameterized as a symmetric TCB function whose param- 
eters are determined from simulated events. The yield of 
'^ISR events in the T{3S) sample, [6.6 ± 1.0(stat.)] x lO'', 
is determined using the off-peak data. A symmetric TCB 
function is used as the PDF for both the xl'^ and ^ ^'^ 



contributions. The peak positions of the Xb compo- 
nents relative to the T'^^^ peak are fixed according to 
the MC simulation. The parameters for the width and 
tail of the TCB function are common to both Xb peaks. 

The K'^ background is modeled using empirical phase 
space functions derived from the MC. Knowledge about 
the ifg transverse momentum distribution is obtained 
from fits to the tt+tt" invariant mass spectrum for the 



Xb 



iv'g-enriched sample, and is used to correct discrepan- 
cies between the data and the MC simulation. The K'^ 
background yield, (348 ± 10) x 10^, is obtained from an 
extrapolation of a fit to the m^j spectrum of the K'^- 
enriched sample, using a scale factor of 2.5 determined 
from MC simulation. The non-peaking background PDF 
is parameterized by a sixth-order polynomial. 

The signal (peaking background) PDF excludes ran- 
dom combinations of tracks that do not originate from 
the signal (background) bottomonium transition. Such 
misreconstructed combinations are included in the non- 
peaking term. 

To improve fit stability, the fit is performed in two 
stages: a preliminary fit to fix background parameters 
followed by a final fit. The peaking background PDF 
parameters and yields are determined from the prelimi- 
nary, x^-based fit in which the signal component is ex- 
cluded from the model. The free parameters in the fit 
are: the yields of the continuum background and the 
peaking background components T^^^, T^~*^, and Xt, 
the continuum background PDF parameters; the overall 
■niR scale of the contribution; the peak positions of 
the T^^"^ and components; the overall widths of 

the PDFs for the T^^'^ , T^^^, and Xb components. 
The x^ per degree of freedom after the preliminary fit is 
364/272 w 1.3, where the largest contributions arise from 
a few isolated bins near 9.79 and 10.02 GeV/c^. As the 
measurement is dominated by systematic uncertainties, 
we evaluate the x^ distribution on simulated pseudoex- 
periments accounting for the dominant sources of system- 
atic uncertainties, and we observe values of x^ greater 
than 364 in more than 7% of the trials. In the final fit, 
all peaking-background parameters except the yields are 
fixed to the values extracted from the preliminary fit. 

The final fit is performed as a scan over the values of 
the hb peak position, with 39 steps in 1 MeV/c^ intervals 
in the range (9.880, 9.920) GeV/c^. At each step, a x^ fit 
is performed for the signal and background yields and the 
continuum background parameters. The fit procedure 
is validated with simulated experiments, and systematic 
uncertainties are evaluated for each point of the scan. 

Figure [T] shows the spectrum and the fit re- 
sult. The non-peaking background component domi- 
nates, with only the prominent T'^^^ and T^^i peaks 
clearly seen above this background. When compar- 
ing the fitted mass of the T'^^^ peak with MC sim- 
ulation, we observe a -H0.44 ± 0.02(stat.) MeV/c^ dis- 
placement in data, which corresponds to a difference 
of 331.50 ± 0.02(stat.) ± 0.13(syst.) MeV/c^ between the 
T(3S') and T[2S) masses. The systematic uncertainty 
is dominated by uncertainties in the lineshape and in 
the track momentum measurement. Details of the lat- 
ter may be found in Ref. 2l|. Other sources of uncer- 



tainty have been investigated and found to be of minor 
significance. These include the fit bias, the cm. boost 
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relative to the laboratory, the background model, and 
the PDF parameters. The position of the T^^^ peak is 
shifted by +1.23 ± 0.02(stat.) MeV/c^ in data with re- 
spect to simulation and corresponds to a difference of 
561.7 ± 0.0(stat.) ± 1.2(syst.) MeV/c^ between the r{2S) 
and T{1S) masses, where the dominant sources of sys- 
tematic uncertainty are the T(2S) momentum in the cm. 
frame and the lineshape model. Figure [2] shows the distri- 
bution of mji after subtraction of the non-peaking back- 
ground. An expanded view of the Xt '"^ region is pre- 
sented in Fig. [31 

The inset of Fig. [2] shows an expanded view of the ht 
signal region. The significance of a signal is evaluated 
at each point of the scan using the ratio, N/cjn, of the 
signal yield N over its uncertainty un- The largest en- 
hancement over background is a 2.2 standard deviation 
(ct) excess (statistical only) at m{hb) « 9.916 GeV/c^. 
Therefore, we do not obtain evidence for an hb signal. 
The fitted h signal yield for m{hb) = 9.900 GeV/c^ is 



-1.1 ± 2.4(stat.)] X 10^ events. Resuhs for the T^" 



product branching fractions are 



T'^\ Xl^\ and xl- 
presented in Table ID 

In the following, reported quantities refer to m{hb) = 
9.900 GeV/c^. The ranges spanned by varying m{hb) in 
the interval [9.880, 9.920] GeV/c^ are given in parenthe- 
ses. The following systematic uncertainties are associated 
with the signal yields. We observe a 10% discrepancy be- 
tween the mji resolution values in data and MC for the 
T^^^ component, and translate this into an uncertainty 
of 0.1 X 10^ (0.0 X 10^ to 0.4 X 10^) events on the K signal 
yield. A systematic uncertainty of 0.4 x 10"^ (0.3 x 10'^ 
to 0.5 X 10^) events is estimated by varying the PDF 
parameters fixed in the fit by ±1 a, varying the overall 



xlO' 




niR [GeV/c1 



FIG. 1: The mR spectrum: data (points) with the fitted 
model (solid line) superimposed. The short-dashed line is the 
contribution from the continuum component. Also shown are 
the r^^^ (dashed curve around 9.79 GeV/c^) and T^^^ (long- 
dashed curve around 10.02 GeV/c^) components. The ht sig- 
nal component is excluded from the superimposed model. 




9.75 9.80 9.85 9.90 9.95 10.00 , 

mR [GeV/c^l 



FIG. 2: The niR spectrum after subtraction of the contin- 
uum background component. The curves represent the fitted 
model (solid), the Tjg^ (dotted), (double-dot-dashed), 
xl'^ (dashed), and xt'^ (dot-dashed) components. Inset: ex- 
panded view in the hb region after subtraction of continuum 
and peaking backgrounds. 
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10.00 



10.01 
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FIG. 3: The ma spectrum in the Xb ' region after subtrac- 
tion of continuum and Ks background components: points 
represent data, while the curves represent the fitted model 
(solid), the x^'^ (dashed), xt'^ (dot-dashed), and Y^^^ (long- 
dashed) components. 



width of the hb PDF by 10%, setting the yield of the 
ISR T{2S) component to ±1 cr of the nominal value, and 
varying the K'^ component normalization and parame- 
ters within their uncertainties. Uncertainties related to 
the continuum background model amount to 0.2 x 10'^ 
(0.0 X 10^ to 0.7 X 10^) events. The additive systematic 
uncertainties on the yields of the T^^^, Xb~^^ and Xb^^ 
components also account for the modeling of the T'^^^ 
tails and for the assumption that the contributions of 
the r{3S) ^ Xxbj'{2P), XbJ'i^P) ^ ^^f:.Ul decay 
chains with J' =^ J or J = are negligible IJ, ll5[ • 

The fit bias on the extracted yields, due to the choice of 
the fit model, is estimated with pseudoexperiments based 
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TABLE I: Summary of results for the signal yields, reconstruction efficiency e, uncertainties on the branching fraction due to 
fit bias, systematic uncertainties on yields and efficiencies, the product branching fraction and the branching fraction for the 
di-pion transition. 





ht [m{ht) = 9.900 GeV/c^] 










Yield 


-1106 ±2432 


31418 ± 1851 


17385 ± 1456 


543839 ± 2928 


906059 ± 7407 




41.8 


25.0 


25.0 


16.7 


47.2 


Fit bias (10~^) 


-0.06 


-0.09 


-0.04 


+0.2 


±0.8 


Yield error (10~^) 


0.01 


0.06 


0.06 


0.6 


0.4 


e error (10~'') 


0.00 


0.05 


0.03 


1.3 


0.8 


ns(io-^) 




1.16 ±0.07 ±0.12 


0.64 ± 0.05 ± 0.08 




17.8 ±0.2 ±1.1 


B (10-3) 


-0.02 ±0.05 ±0.06 


9.2 ±0.6 ±0.9 


4.9 ±0.4 ±0.6 


30.0 ±0.2 ±1.4 





on fully SHTLulated Monte Carlo samples. We estnTiate a 
fit bias on the h signal yield of -2.8 x lO'"^ (-3.0 x 10^ 
to ±0.4 X 10^) events. Fit biases for the other di-pion 
transitions are listed in Table H] We do not correct the 
signal yields but rather assign the bias as a systematic 
uncertainty. 

The following systematic uncertainties are associated 
with the reconstruction efficiency e. The uncertainty due 
to the track- reconstruction efficiency is 3%. To assess the 
impact of data-MC differences on the tt+tt" candidate se- 
lection efhciencies, we compare the relative variations of 
the T^^^ yield in data and MC when excluding selection 
requirements one at a time, and assign the full observed 
discrepancy to the systematic uncertainty. A total un- 
certainty of 2.3% in e is obtained, including the statis- 
tical uncertainty (0.6%) in the T^~^'^ yield. The uncer- 
tainty in the number of T{3S) events is 1.1%. The above 
multiplicative systematic uncertainties affect the prod- 
uct branching fractions of all di-pion transitions studied 
in this analysis. Differences in the selection efficiencies 
resulting from different angular distributions of the hb 
decay products and different hb hadronization models in 
the MC simulation contribute a 5% uncertainty. Model 
uncertainties in the simulation of the di-pion kinemat- 
ics, bottomonium hadronization, and T(2S) production 
channel (where applicable) in the T^-*\ j'3-i-2 p2]^ xl'^^^ 
and Xb~^^ decay chains result in systematic uncertainties 
on the efficiency of 1.3%, 0.5%, 0.6%, and 0.6%, respec- 
tively. 

Product branching fractions are calculated by divid- 
ing the fitted yield by the efficiency and the number 
of T{3S) events, and are summarized in Table HI For 
m{hb) = 9.900 GeV/c^ we find the branching fraction 
Br = B[T{3S) Tr+Tr^hb] = (-0.2 ± 0.5 ± 0.6) x IQ-^, 
where the first uncertainty is statistical and the second 
systematic, and set an upper limit (UL) Br < 1.2 x lO^'* 
at 90% CL. The UL is calculated assuming a Gaus- 
sian sampling distribution f{Br) for Br, which accounts 
for statistical and systematic uncertainties, and deter- 
mining the value of UL for which Jq'" f{Br)dBr ~ 
0.9 X Jg°° f{Br)dBr- Figure [3] reports the branching 




9.880 9.885 9.890 9.895 9.900 9.905 9.910 9.915 9.920 

m(h ) [GeV/c'] 

b 



FIG. 4: Results for the T{3S) — -K^n^hb branching frac- 
tion (points with statistical and systematic uncertainties) as 
a function of the assumed hi, mass. The triangles indicate the 
upper limits at 90% CL. 



fractions Br and the 90% CL ULs as a function of 
the assumed hb mass. The branching fractions of the 
X6i(2P) -± 7r+7r-X6i and Xfc2(2P) -J> 7r+7r-X62 tran- 
sitions, given in Table U are derived by correcting for 
the branching fractions of the T{3S) 7Xbi(2P) and 
r(35) 7X62 (2P) decays 4, 21, respectively. The ex- 
tracted values are in reasonable agreement with those 
found in the study by the CLEO Collaboration 0, [3] , 
where the two transitions could not be separated exper- 
imentally. 

In summary, we present an inclusive analysis of the 
TT+TT^ recoil mass spectrum in T(3S) decays. We measure 
the branching fraction 

B[r{3S) Tr+TT-T{2S)] 

= (3.00 ± 0.02(stat.) ± 0.14(syst.))%. 

This value is in reasonable agreement with, and more 
precise than, the current world average (2.45±0.23)% 

The measured T{3S)-T{2S) mass difference is 331. 50± 
0.02(stat.) ± 0.13(syst.) MeV/c^. 
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We extract the product branching fractions 

B[Ti3S) -> Xxbi{2P)] X B[xbi{2P) ^ Tr+n-Xbi] 
= (1.16±0.07±0.12) X 10"^ 

B[r{3S) ^ Xxb2{2P)] X B[Xb2{2P) ^ 7T+TT-Xb2] 

= (0.64 ±0.05 ±0.08) X 10"^ and 

B[r{3S) Xr{2S)] X B[r{2S) n+n-r] 
= (1.78 ±0.02 ±0.11)%. 

A search for the state, the ^Pi state of bottomo- 
nium, in T(3S') Tr^w^hf, decays does not provide 
evidence for this decay mode, and assuming the hb 
mass to be 9.900 GeV/c^, we set a 90% CL upper limit 
Br < 1.2 X 10-"*. We exclude, at 90% CL, values of 
Br above 1.8 x 10^^ for a wide range of assumed hb 
These results disfavor the calculations of 
Similarly, a recent measurement of the 



mass values 
Refs. [i-ini 

T{l^Dj) T{lS)TT+n- branching fraction |2| disfa 



inc hing 

vors the calculations of Ref. 13, Hf • The predictions of 
Ref. are at least one order of magnitude smaller and 
are not contradicted by our result. 
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